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C/3 ■ ABSTRACT 
a _ 

We present analysis of both the resolved X-ray emission line profiles and the 
broadband X-ray spectrum of the 02 If* star HD 93I29A, measured with the Chandra 
^ [ HETGS. This star is among the earliest and most massive stars in the Galaxy, and 

\^ ■ provides a test of the embedded wind shock scenario in a very dense and powerful wind. 

OO ' A major new result is that continuum absorption by the dense wind is the primary 

. cause of the hardness of the observed X-ray spectrum, while intrinsically hard emission 

' from colliding wind shocks contributes less than 10% of the X-ray flux. We find results 

I consistent with the predictions of numerical simulations of the line-driving instability, 

, including line broadening indicating an onset radius of X-ray emission of several tenths 

i?*. Helium- like forbidden-to-intercombination line ratios are consistent with this onset 
radius, and inconsistent with being formed in a wind-collision interface with the star's 
closest visual companion at a distance of 100 AU. The broadband X-ray spectrum is 
fit with a dominant emission temperature of just kT = 0.6 keV along with significant 
, wind absorption. The broadband wind absorption and the line profiles provide two 

^ ' independent measurements of the wind mass-loss rate: M = 5.2^}'5 x 10~^ Mq yr~^ 

■ ■ ■ and M — 6.8^2 2 ^ 10^^ Mq yr^^, respectively. This is the first consistent modeling 

of the X-ray line profile shapes and broadband X-ray spectral energy distribution in a 
massive star, and represents a reduction of a factor of 3 to 4 compared to the standard 
Ha mass-loss rate that assumes a smooth wind. 

Key words: stars: early-type - stars: mass-loss - stars: winds, outflows - stars: 
individual: HD 93129A - X-rays: stars 



1 INTRODUCTION 

With a spectral type of 02 If (|Walborn at al.ll2002l ). HD 
93129A is among the earliest, hottest, most massive and 
luminous O stars in the Galaxy. As such, it has an ex- 
tremely powerful wind, with a terminal velocity in excess 
of 3000 km s~^, and a m ass-loss rate thought to be in 
excess of 10'^ Mq yr~^ (|Repolust et all |2004 ). In fact, 

* E-mail: colien@astro.swarthmore.edu 



HD 93129A has been considered to ha ve the highest mass- 
loss rate of any O star in t he Galaxy l|Taresch et al.|[l997l : 
iBenaglia fc Koribalskil [20041 ) . This star, therefore, provides 
an interesting test of the wind-shock paradigm of massive 
star X-ray emission, both because of the tremendous kinetic 
power in its wind and because of the dominant role played 
by X-ray absorption in such a dense wind. In order to study 
the X-ray emission and absorption we have analyzed the 
Chandra HETGS spectrum with a focus on individual line 
profiles. In addition, we analyze the low-resolution, zeroth- 
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order ACIS CCD spectrum in order to determine the rel- 
ative contributions of high-temperature thermal emission 
and wind attenuation to the observed spectral hardness. 
These complementary analysis techniques provide informa- 
tion about the temperature and kinematics of the X-ray 
emitting plasma, about its spatial distribution, and about 
the wind mass-loss rate. 

X-ray emission from O stars is attributed to three 
mechanisms: (1) Embedded Wind Shocks (EWS), gener- 
ally assumed to be associated with t he Line-Driving In- 
stability (LDI ) ("Lucy & White" 'igSO"; 'Owoc ki et all Il988l : 
iFeldmeier eTal. 1997 ; Kahn_et_al, 2001) ; (2) Colliding Wind 
Shocks (CWS) in some binary svstems Jstevens et al.l 1 19921 : 
lAntokhin et al. l l2004l : lPittard fc Parkinlf2010l ); and (3) Mag- 
netically Confined Wind Shocks (MCWS) for stars with sig- 
nificant dipole magnetic fields ("Babel fc Mon tmerld 1 19971 : 
lud-Doula fc Owocki|[2003 ; [Cagnc ct al. 2005). Of these, the 
EWS mechanism is assumed to operate in all O stars, while 
CWS may dominate in massive binaries with strong enough 
winds and MCWS in those small number of O stars with 
strong, large-scale magnetic fields. The EWS mechanism 
produces plasma of several million degrees and is associ- 
ated with relatively soft X-ray emission, while the other two 
mechanisms produce stronger shocks, higher temperatures, 
and harder X-ray emission. However, it should be kept in 
mind that soft X-ray absorption by the bulk wind can harden 
the observ ed X-rays from EWS in O stars with high mass- 
loss rates l|Leutenegger et al]|2010l ). 

In high-resolution X-ray spectra, the hall mark of em- 
bedded wind shocks is bro ad emission lines (|Kahn et al.l 
I2OOII : ICassinelh et al.ll200ll ). By analyzing the widths and 
profile shapes of individual X-ray emission lines in the grat- 
ing spectra of O stars, the kinematics of the hot. X-ray emit- 
ting plasma embedded in the warm, partially ionized bulk 
wind can be determined, testing the predictions of the EWS 
scenario. Furthermore, due to preferential absorption of red- 
shifted line photons from the far hemisphere of O star winds. 
X-ray emission lines from embedded wind shocks have a 
characteristic blue-shifted and skewed sh ape, in proportion 
to th e characteristic wind optical depth (jOwocki fc CohenI 
[2003). It recently has been shown that a wind mass- loss rate 
can be determined by fitting the ensemble of derived charac- 
teristic o ptical depths, given a model of the bulk wind X-ray 
opacity (|Cohen et al.l 120101 ). The initial application of this 
technique to Pup (04 If) provided a mass- loss rate deter- 
mination that was roughly a factor of three lower than the 
traditional value derived from the strength of the Ha emis- 
sion under the assumption of a smooth, unclumped wind. 
This l ower value is co nsistent with other recent reassess- 
ments (|Puls et al.ll2006l ) that do account for small-scale wind 
clumping, which affects density-squared diagnostics such as 
Ha emission strength. The associated clumping factors are 
consis tent with those seen in numerical simulations of th e 
LDI (jRunacres fc Owockil I2OO2I : iDessart fc Owockil l2005h . 
Indeed, there is a consensus emerging that many, if not all, 
O stars' mass- loss rates must b e lowered by factors o f several 
due to the effects of clumping l|Hamann et al ][200i). In this 
context, the X-ray line profile mass- loss rate diagnostic is 
especially useful, as it is not affected by small-scale clump- 
ing, as long as the individ ual clumps are optically thin to 
X-rays (|Cohen et al.ll2010l ). 

Both numerical simulations and the lack of significant 



observed X-ray variability indicate that clumps in O star 
winds are on quite small scales, with sizes ^ <C -R*. Since even 
the entire wind of HD 93 129 A is only marginally optically 
thick to bound-free absorption of X-rays, it is very likely that 
such small clumps will be individually quite optically thin 
to X-rays. This means they cannot have much of the self- 
shadowing that would reduce the exposure of wind material 
to X-rays, and would thus lead to a significant "porosity" 
reduction in the overall absorption. An important point here 
then is that, while such porosity requires a strong clumping 
(with individual clumps that are optically thick) , a clumped 
wind need not be porous (if the clumps are optically thin)Q 
HD 93129A lies in the Trumpler 14 cluster in the Carina 
nebula, at a distance of 2.3 kpc, but with a modest visual 
extin ction of Ay — 2.3 (|Townslev et al.ll201ll : ICagne et al.l 
I2OIII ). It has a visual companion, HD 93129B, at a sep- 
aration of 2.7", with a spectral type of 03.5, and it also 
has a closer visual companion (HD 93129Ab) detected us- 
ing the H ST Fine Guidance Sensor (FGS) at a separation 
of 0.053" (|Nelan et al.ll2004 I20T0I ). This closer companion 
is also estimated to have a spectral type of 03.5. Non- 
thermal radio emission has been detected from the system, 
presumably indicatin g the existence of colliding wind shocks 
(|Benaglia et al.ll2006l ). This motivated the initial CWS in- 
terpretation of low-resolution Chandra CCD spectral mea- 
surc ments of the r elatively hard X-rays from HD 93 129 A 
( Evans et al.ll2003l ). although given the large separation of 
components Aa and Ab (over 100 AU « 1000 R*), any X-ray 
emission associated with the wind interaction zone should be 
relatively weak and the overall emission is likely dominated 
by EWS emissi on arising mu c h clos er to the photosphere of 
component Aa. ICagne et al.l (|201ll ) has recently suggested 
that only 10 to 15 percent of the X-ray emission is due to 
CWS X-rays. Furthermore, given the high wind mass-loss 
rate of HD 93129A, significant attenuation of the soft X-ray 
emission is expected. This paper represents an attempt to 
account for wind absorption when analyzing the star's X-ray 
properties. 

Detailed analysis and modeling of the o ptical and UV 
spectr a of HD 93 129 A was presented by iTaresch et al.l 
([1993). That work included an analysis of the interstellar 
absorption features in the UV, which yielded a hydrogen 
column density measurement of TVh = 2.5 x 10^^ cm~^ that 
is basically consistent with the observed reddening and in- 
ferred extinction (A^h = 3.7 x 10^^ cm~^ is implied by the 
color excess). Important results from the UV and optical 
spectral analysis include a very high wind mass-loss rate 
and terminal velocity in addition to a bolometric luminosity 
in excess of 2 x 10^ Lq and evidence for non-solar abun- 
dances in line with CNO processing (a nitrogen abundance 
several times solar and carbon and oxygen abundances a 

^ In this paper we refer to "optically thin clumping" when dis- 
cussing clumps that do not have a porosity effect but do have 
an effect on Ha, reserving "porosity" to describe the effects of 
optically thick clumps. For the unmodified term "clumping" , the 
reader should bear in mind that density-squared diagnostics such 
as Ho will still always be affected, but that X-ray transmission 
will be affected by an associated porosity if and only if the clumps 
are optically thick. Finally, we note that a given clump or clump 
distribution may be optically thin at one wavelength and thick at 
others. 
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Table 1. S tellar and wind pa rameters adopted from 
iTaresch et"aD l(l997i ) and lRepolust et al. (,2001 ) (in parenthe- 
ses). 



parameter value 



Mass (95 Mq) 

Teff 52000 K (42500 K) 

19.7 R0 (22.5 R0) 
log L, 6.4 L0 (6.17 Lq) 

Drotsini (130 km s~^) 

3200 km s'^ (3200 km s"!) 
/3 0.7 (0.8) 

M 1.8 X 10-5 Mq yr-i (2.6 X 10"^ Mq yr"!) 



factor of a few lower than solar). These authors estimate a 
zero-age main sequence mass of 120 M f ^^ for HP 93129A. 
More recent analysis by iRepolust et al.l ()2004l ) finds a sig- 
nificantly lower effective temperature and a modestly lower 
bolometric luminosity along with a higher wind mass-loss 
rate. The stellar parameters derived from these two studies 
are summarized in Table[T] We note that the presence of the 
close binary companion will affect the radius determination, 
and related quantities. Given that the companion is about 
a mag nitude dimmer than the primary, this effect will be 
smaU (|Repolust et al.ll2004l ). 

We describe the Chandra HETGS data in §2. In §3 we 
analyze the line profiles seen in the grating spectrum and 
also fit global thermal emission models to the low-resolution 
but higher signal-to-noise zeroth order spectrum. We discuss 
the implications of these analyses for the three mechanisms 
of O star X-ray emission in §4, and summarize our conclu- 
sions in 55. 



2 THE Chandra DATA 

The data we use in this paper were taken between 8 Novem- 
ber 2005 and 5 December 2005 in seven separate point- 
ings, with a total effective exposure time of 137.7 ks. All 
observations employed the Advanced CCD Imaging Spec- 
trometer with the High Ener gy Transmission Grat ing Spec- 
trometer (ACIS-S/HETGS) (|Canizares et al.l l2005'). provid- 
ing dispersed spectra in the MEG and HEG grating arrays, 
as well as a low-resolution CCD spectrum from the zeroth- 
order i mage. The grating data have been pres e nted previ- 
o uslv in IWestbrook et all ((200^) and IWalbornI (|2008l ). and 
in lNichols et al.l (120111 ^ 

The field is crowded with X-ray sources, but only HD 
93129B, our target's visual companion, is bright enough and 
close enough to HD 93129A to pose a potential problem in 
the data extraction. In Fig.[T]we show the center of the ACIS 
detector for the longest of the seven separate exposures, with 
the MEG and HEG extraction regions indicated, and the 
zeroth-order images of components A and B labeled. The 
CCDs that compose the ACIS detector have modest intrin- 
sic energy resolution (of roughly E / AE — 20 to 50), allowing 
us to assess the spectral energy distribution of each source 
as indicated by the color coding in the figure. For this par- 
ticular exposure, component B, which is 2.7" southeast of 
component A, lies in the MEG and HEG extraction regions. 



Figure 1. The central region of the ACIS detector from the 
longest single exposure (Obs ID 7204; exposure time of 34 ks, 
corresponding to 25 percent of the total exposure time), showing 
the zoroth order image of HD 93129A (indicated by the circle), 
along with several other sources. HD 93129B is labeled and lo- 
cated 2.7" to the southeast of component A (and indicated by 
the modestly elongated ellipse). The detected photons are color 
coded according to energy, with low energies (0.5 to 1.5 keV) red, 
medium energies (1.5 to 2.5 keV) green, and high energies (2.5 
to 8 keV) blue. Component B is clearly separable from compo- 
nent A and is softer and weaker (by roughly a factor of ten) than 
component A. Its relative weakness and softness are quantified 
in Fig. [2] The negative and positive MEG and HEG grating arm 
extraction regions are outlined in white (HEG is more vertical, 
MEG is angled to the upper right and lower left). Note that these 
have a somewhat different orientation in each observation. The 
zeroth order image of component B lies mostly within both of 
these extraction regions, indicating possible contamination of the 
dispersed spectra of HD 93129A. 



and therefore is a potential source of contamination of the 
dispersed spectra of component A. 

Due to the varying roll angle of the Chandra instrument, 
component B does not lie fully in the extraction regions of 
component A for all of the seven separate observations. Fur- 
thermore, as shown in Fig. [2] its photon flux is an order 
of magnitude lower than that of component A. Finally, the 
spectral energy distribution of component B is much softer 
than that of component A (this can also be seen, qualita- 
tively, in the color coding of Fig. [TJ. In fact, there are al- 
most no counts from component B at energies greater than 
1.2 keV, corresponding to a wavelength of 10 A. Thus, the 
dispersed spectra of component A are unaffected by con- 
tamination from component B below this wavelength. At 
wavelengths above 10 A, however, there is likely to be mod- 
est contamination. We note that the 2.7" offset corresponds 
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Energy (keV) 



Figure 2. The extracted ACIS CCD spectra — zeroth order - 
from the seven coadded pointings (circles) and the same for HD 
93129B (squares). Note that above 1.2 keV (below 10 A), the 
contamination of the grating spectra of component A by that of 
component B should be negligible. 



Mg xil), especially for silicon. This qualitative impression 
is borne out by quantitative modeling presented in §3.3. 

Only a handful of lines are present in the MEG spec- 
trum (and even fewer in the lower signal-to-noise HEG spec- 
trum). The small number of visible lines is affected by the 
overall low signal-to-noise and the aforementioned attenua- 
tion, which renders the normally quite strong Fe xvil and O 
lines longward of 15 A completely absent. And it is exacer- 
bated by the very large line widths, which spread the modest 
number of line photons over many pixels. After discarding 
the very weak lines longward of 10 A due to contamination 
from component B, we are left with five lines and blended 
line complexes detected with greater than 3a" significance: 
the He-like S xv complex near 5.1 A, the Lya line of Si xiv 
at 6.18 A, the He-like Si xill complex near 6.7 A, the Mg 
xil Lya line at 8.42 A, and the He-like Mg xi complex near 
9.2 A. The weak Mg xi He/3 line at 7.85 A is not detected 
above the 3a threshold. Not shown in the figure are very 
weak Ne ix and Ne x lines between 12 A and 14 A, which 
are affected by contamination by the softer X-ray emission 
from component B, as discussed above. 



to a relative shift in the dispersed MEG spectrum of roughly 
1500 km s~^, which could artificially broaden the emission 
lines of component A. We therefore restrict the analysis of 
the grating spectra in this paper to wavelengths shortward 
of 10 A. 

We note also that the close binary companion Ab de- 
tected at 0.05"with the FGS is completely unresolved by 
Chandra so the spectra we analyze in this paper are a com- 
posite of the two components. However, the spectral type 
of Ab is tentatively ident ified as b eing the same as that of 
component B (03.5 V) ((Nelan et al.i i2004i). so the compos- 
ite spectrum, especially below 10 A, is likely dominated by 
component Aa. We discuss the colliding wind X-ray emission 
due to this close visual companion in §3.3. 

After centroiding the zeroth-order image of component 
A in each of the seven separate observations, we coadded the 
observations and extracted the first-order MEG and HEG 
spectra, as well as the zeroth-order ACIS CCD spectrum 
shown in Fig. [S] The coadded negative and positive first- 
order spectra (both MEG and HEG) are shown in Fig. [S] 
The MEG spectrum has a FWHM resolution of 2.3 mA and 
the HEG, with a lower sensitivity, has a resolution of 1.2 
mA. In the 5 to 10 A region where most of the counts are, 
these correspond to resolving powers of roughly A/ A A — 300 
to 600 in the MEG. 

Even with the coaddition of seven separate pointings, 
for an effective exposure time of 137.7 ks, the spectra have 
quite low signal-to-noise. This is due both to the great dis- 
tance to HD 93129 (2.3 kpc) and also to the large interstellar 
column density, which causes significant attenuation of the 
soft X-ray emission from the star. As we will show in !j3l 
there is also significant X-ray attenuation from the dense 
stellar wind. The appearance of the Chandra grating spec- 
tra of HD 93129A is certainly significantly harder than that 
of other normal O stars, but this is due primarily to atten- 
uation and not to high temperatures, as can be seen qual- 
itatively in Fig. [3] from the dominance of the He-like lines 
(Si xill, Mg xi) over the corresponding H-like lines (Si xiv, 



3 SPECTRAL ANALYSIS 
3.1 Resolved emission lines 

By fitting a simple , empirical line profile model 
(|Owocki fc CohenlboOll ) to the Doppler-broadened emission 
lines, we can simultaneously determine the kinematics of 
the X-ray emitting plasma and the degree of attenuation 
by the wind in which the hot, shock- heated plasma is 
embedded. The specific parameters of the lOwocki fc CohenI 
(|200lf ) model are the onset radius of the X-ray emission 
(Ro) and the fiducial optical depth of the bulk wind, 
r, = KM/47rR»Voo. Note that r, is expected to vary from 
line to line due to the wavelength dependence of the wind 
opacity, k. 

The fitted -Ro values derived from the data are ex- 
pected to be several tenths of a stellar radius above the 
photosphere, ba sed on simulations of the line-driving in- 
stability (LDI) (lOwocki et al.l Il988l : iFeldmeier erakl Il997l : 
iRunacres fc Qwockill2002l) ." By deriving values for this pa- 
rameter from the individual lines in the Chandra spectrum 
of HD 93129A, we can test the LDI scenario for embedded 
wind shocks in the most extreme O star wind. And the val- 
ues we derive for the fiducial optical depth, r, , can be used 
to derive a mass-loss rate by fitting the ensemble of values, 
given a model of t he wind opacity, as has been shown by 
ICohen et al.l (|2010D . 

F ollowing the procedure described in I Cohen et al.l 
(|201(]t ). we assess the continuum level near each line by fit- 
ting a small region of the spectrum on either side of the 
line, and then fit a profile model plus the continuum model 
(with the level fixed at the value found from fitting the 
nearby continuum). We allow the normalization, Ro, and 
r«, to be free parameters of t he fit, while fixing t he veloc- 
ity law parameter at = 0.7 (|Taresch et alj[l997l) and the 
terminal velocity at the value determined fr om the analy- 
sis of UV observations, Vo o = 3200 km s~^ (|Taresch et al.l 
ll997l : lRepolust et al.|[20o3 ). We find t he best-fit m odel pa- 
rameters by minimizing the C statistic (|Cashll 19791 ). and as- 
sign confidence limits individually to each model parameter 
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Figure 3. The extracted MEG (top) and HEG (bottom) spectra from the seven coadded pointings. Note the different y-axis scales on 
the two figures. The wavelengths of lines expected to be present in normal O star Chandra spectra are indicated by the vertical dotted 
lines. 



(while allowing the other free paramet ers to vary) according 
to the formahsm in lPress et al] (|2007l '). We fit the MEG and 
HEG data simultaneously. We perform all of this modeling 
and data analysis in xspec v. 12.6, using the custom model 
windprofil^ 

For the helium-like complexes, which effectively com- 
prise three blended lines each (resonance (r), intercombi- 
nation (i), and forbidden (f)), we fit three profile models si- 
multaneously (using the custom model hewind), with the Ro 
and r* parameters tied together for each of the three lines, as 
described in Lcutcnc g gcr c t al. (200g). The overall normal- 
ization and the Q = {f+i)/r ratio are explicit fit parameters, 
while the diagnostically important TZ = f /i ratio is a func- 
tion of radius, via the radial dependence of the photoexcita- 
tion rate of electrons out of the upper level of the forbidden 
line into the upper level of the intercombination line. This 
physics is controlled in our modeling by the same Ro param- 
eter that describes the onset radius of X-ray emission. Lower 
values of -Ro give more plasma close to the photosphere, 
where it is more strongly affected by photoexcitation, pro- 
viding a lower TZ = f /i value than if Ro were larger and the 
overall photoexcitation rate were lower. Thus, Ro controls 
both the line widths and the relative strengths of the inter- 
combination and forbidden lines in the observed spectrum 
in a self-consistent manner. For the photoexcitation model- 
ing, we use UV fluxes from a TLUSTY Teff = 425 00 K, log 
g = 3.75 model atmosph ere (|Lanz fc Hubenvll200^ ). and the 
atomic parameters from lDere et al. I (|2007h . 

The number of free parameters in the models - both for 
the single lines and for the He-like complexes - is kept to 



The windprofile model's implemen- 

tation in XSPEC is described at 

heasarc.gsfc.nasa.gov / docs / xanadu/xspec / models /windprof .html, 
as is the hewind model we use to fit helium-like complexes. 



a minimum. We fix the wind terminal velocity {voo ~ 3200 
km s~^) and the velocity law parameter (/3 = 0.7) at the 
value determined from the UV data. There is an extensive 
discussion of the sensitivity of the important model param- 
eters, r« and Ro, to these, and other, fixed parameters as 
well as the effect of such factors as background subtrac tion 
and continuum placement in §4.3 of ICohen etHI l|2010l ), to 
which we refer the reader. 

The results of our wind profile fits to the five lines and 
line complexes are shown in Fig. |4] for the single Lya lines 
and Fig.[S]for the three He-like line complexes. The quantita- 
tive results are summarized in Tab. [2] We note that when we 
increase the wind velocity parameter, /3, from 0.7 to 1.0, the 
characteristic optical depth values, r*, increase by roughly 
30%, while the onset radii, Ro, increase by several tenths of a 
stellar radius. Similarly, there is sensitivity of the important 
derived parameters to the assumed wind terminal velocity. 
For a higher terminal velocity of 3400 km s~^, the best- 
fit r* values decrease by roughly 20% and increase by the 
same amount when we use a lower terminal velocity of 3000 
km s~^. The onset radii, Ro, vary by roughly 10% for these 
changes in terminal velocity. We note that these systematic 
uncertainties - especially for r, - are small compared to the 
statistical errors. 

Before we move on to interpreting these model-fitting 
results, we report on a few experiments involving fitting dif- 
ferent types of models to these same emission lines. Specif- 
ically, Gaussian profiles have traditionally been fit to the 
broadened emission lines seen in O stars. Therefore, we fit 
the highest signal-to-noise single line in the spectrum, the 
Mg xil Lya line, with an unshifted Gaussian. As expected 
from the asymmetry seen in the bottom panel of Fig. [4] the 
fit is poor (the wind profile fit is preferred at > 99.99%). Al- 
lowing the Gaussian centroid to be a free parameter, we 
find an improved fit with a large centroid blue shift (to 
A = 8.401 A), equivalent to -730 ± 200 km s"\ and a 
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Table 2. Wind profile model fit results 



ion wavelength*^ t, Ro normalization'' 

(A) ~ (R*) (10-^ ph cm-2 r-l) 

S XV 5.0387, 5.0648, 5.1015 0.32+^3^*^ l-l^^ i* 4.11+j 

Si XIV 6.1822 l-64ti:26 l-Wt'o? l-48l:^^ 

Si XIII 6.6479, 6.6866, 6.7403 0.8lt || 1.37^ 16.8tJ:g 

Mg XII 8.4210 l-OSt^gg^ l-44l;27 'i-OOt'l^ 

Mg XI 9.1687, 9.2297, 9.3143 2.84l;^''^g l-SOt ^g 16. 1^^;^ 



Closely spaced doublets in the Lyo lines and He-like intercombina- 
tion lines are fit with a single profile model centered at the emissivity- 
weighted wavelength of the two components. 

For the He-like complexes, the total normalization of all the lines 
in the complex is indicated. 




6.15 6.20 
Wavelength (A) 



6.25 




8.35 



8.40 _ 8.45 

Wavelength (A) 



8.50 



Figure 4. The MEG data for both individual lines with sufficient 
signal-to-noise to warrant profile fitting, along with the best-fit 
profile model for each line (histogram). These are the Lya lines 
of Si XIV, at 6.18 A (top), and Mg XII, at 8.42 A (bottom). The 
laboratory rest wavelength of each line is indicated by a vertical 
dashed line and the Doppler shifts associated with the (positive 
and negative) terminal velocity are indicated by the vertical dot- 
ted lines. Poisson error bars are indicated on each data point. Note 
that the best-fit model for the Mg XII line is based on jointly fit- 
ting the HEG and MEG data, though we show only the MEG 
data here. 



width of o- = 1200li?^ km s"^ This shifted-Gaussian fit 
is statistically indistinguisfiable from the wind profile model 
fit. Although the wmdprofile fit is more meaningful, as it is 
based on a physically realistic model while having no more 
free parameters than the Gaussian model, the Gaussian fit 
confirms and quantifies the large line widths and blue shifts 
expected from EWS emission. 

Along similar lines. He-like complexes can be fit with 
three Gaussians, providing a direct measure of the f/i ra- 
tio, from which a single radius of formation can be inferred. 
We show these results in Fig. |6l For the Mg xi complex we 
find a low ratio of 0.40l'^y. We next fit the same model, 
but with the f/i ratio fixed at the "low density limit" that 
would be expected if the X-ray plasma were very far from 
the photosphere, where photoexcitation cannot alter the ra- 
tio. This fit (which is also shown in Fig. O is poor compared 
to the f/i — 0.4 fit, implying that the hot plasma is rela- 
tively close to the photosphere. This conclusion is in good 
agreement with the result from fitting the three wind-profile 
model (hewind), in which the onset radius of a distributed 
source of X-ray emitting plasma, Ro, is constrained to be 
within 1 stellar radius of the photosphere. 

Porosity, due to optically thick clumps, can affect X- 
ray li ne profiles by reducing the average optical depth of the 
wind (jOskinova et al.ll2006l ). and is another effect we might 
consider when fitting the emission line data. Ifowever, the 
porosity length that is necessary to provide a measurable 
effect is quite large (|Owocki fc CohenI |2006| ) compared to 
the small-scale structure in state-of-the-art 2-D radiatio n 
hydrodynamics simulations (|Dessart fc Owockill2003l , |2005| ) . 
To explore the effects of porosity on the fit quality and on 
the other parameters, we fit the Mg xil Lya line with a 
model having an effective opacity modified by porosity from 
sp herical clumps. This mo del is similar to the one described 
in lOwocki fc CohenI l|2006l ) - using the same porosity- length 
{h{r) = £/ f, where £ is the characteristic clump size and / 
is the clump filling factor) formalism - but here employing 
a radial clump distribut ion determined by th e wind beta- 
velocity law, as used by :Oskin ova et al] l|2006l) . We fix the 
characteristic optical depth at r» = 4.2, which is the value we 
would expect at the wavelength of the Mg L ya line assum- 
ing a mass-loss rate of 1.8 x 10~^ Mq yr"'^ (|Taresch et al.l 
[l993). The fit we obtain is similar in quality to the windpro- 
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Figure 5. The MEG data for the hehum-hke hne complexes: S 
XV (top), Si XIII (middle) and Mg XI (bottom). We show line 
center and terminal velocity indicators as in the previous figure, 
but here, we show three groups; one for each of the three lines 
in each complex. The best-fit models are based on jointly fitting 
the HEG and MEG data (except for the weaker S XV complex, 
where we fit only the MEG data). Note the very weak forbidden 
line in both the Mg and Si complexes (longest wavelength of the 
three lines in each complex). 
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Figure 6. The MEG data for the Mg XI line complex fit with 
three shifted Gaussians. When the f /i ratio is a free parameter 
(red), it is low (0.4^'^i^). When the ratio is fixed at the normal, 
equilibrium limit (//i = 2.7) with no photoexcitation, as would 
be expected far from the star (green), the fit is poor. 



file fit reported in Table [2] But in order to achieve this good 
fit, a very large terminal porosity length of ftoo = 5.3 i?* 
is required (with a 68% lower confidence limit of hoo = 2.5 
Rt), where the radially varying porosity length is given by 
h{r) = /ioo(l - R*lrf. If we use 2.6 x 10 "'^ yr~^ as 
the s tandard, smooth-wind mass-loss rate l|Repolust et al.l 
12004 ). then the required porosity length is even larger. Even 
the minimum porosity length of hoo = 2.5 is inconsistent 
wit h the numerical simulations of the line-driving instabil- 
ity (|Dessart fc OwockilbOoA |2005| ). requiring, for example, 
clumps that are individually 0.25 i?« in scale in a wind with 
a uniform filling factor of /v = 0.1. 



3.2 Mass- loss rate determination from the line 
profiles 

Given the characteristic optical depth values (r«, listed in 
Table [2ll obtained from fitting each line with the windprofile 
model, we can derive constraints on the wind mass-loss rate 
of HD 93 129 A. The characteristic optical depth is defined as 
r, = KM/47rR,Voo, so that by assuming a model of the bulk 
wind's X-ray opacity, k(A), as well as a value for the stellar 
radius and wind terminal velocity, we can fit the wavelength- 
dependent r* values with the mass-loss rate as the only free 
parameter of the fit, as M = 47rR*Voort(A)/K(A). We have 
demonstrated this process for analyzing the Cfeandra gratiri g 
spectrum of the O superg iant C Pup (|Cohen et al.l l201Gl ). 
Because of the lower quality of the Chandra spectrum of HD 
93129A compared to that of Pup, and because of the much 
smaller number (5 vs. 16) and narrower wavelength range 
of usable lines, the wavelength trend in r* is not apparent 
for HD 93129A. However, the five r* values are certainly 
consistent with the expected wavelength trend. 

We constructed a model of the bulk wind opacity for 
HD 93129A, assuming that H and He are fully ionized, and 
assuming typical values for the ionization balance of metals 
(generally dominated by triply io nized states). Furtherm ore, 
we assume solar abundances from lAsplund et all (j2009h . ex- 
cept for C, N, and O, which are altered by CNO processing. 
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Figure 7. The five r* values fit with a single opacity model in or- 
der to derive M. The best-fit r* model (M = 6.8 X ID"® Mq yr-^) 
is shown as a solid line, while the dotted line represents that form 
of r* that would be e xpected if the traditi onal mass-loss rate of 
2.6 X 10"^ M0 yr-l l lRepolust et al.ll2004^ were correct. 



according to the spectral analysis of iTaresch et al.l l|l997l ). 
We take the C, N, and O values from these authors, but 
rescale them so that the sum of the abundances of these 
thr ee elements is e qual to the sum of the C, N, and O in 
the lAsplund et af] (120091 ') solar abundances. This effectively 
gives us Zc = 0.25 , Zn = 3.02, and Zq = 0.47, where 
I Asplund et al] l|2009l ) is the solar reference. We should note, 
however, that because the emission lines we derive r* val- 
ues from are all at short wavelengths, the alterations to C, 
N, and O abundances have almost no effect on our results. 
Changing the overall metallicity - which we assume to be 
solar - would have an effect, however. The derived mass-loss 
rate scales inversely with the metallicity. 

The result of fitting the r, values for the mass-loss rate 
is shown in Fig. [T] The best-fit mass-loss r ate is 6.8t.o'^ x 
10~^ M yr~^, using the stellar radius from Ropolust ct al.' 
l|2004h . This represents a fact or of three ( Tarcsc h et al- f 997i) 
or four l|Repolust et al.|[2004l ') reduction with respect to tra- 
ditional Ha-based determinations that assume a smooth 
wind (and thus ignore optically thin clumping). 

To check for consistency with the observed Ha, we 
modeled the profile with the line-blanketed, non-LTE, uni- 
fied (photosphere -l-wind) model atmosphere code fastwind 
l|Puls et al.ll2005l '). which accounts for optically thin clump- 
ing (using the filling factor approach) in the calculations 
of the occupation numbers and the correspo nding synthetic 
spectr a. Stellar parameters were taken from lRepolust et al.l 
l|2004 ). where the observational data also are described. 
However, we adopt the mass-loss rate derived from our X-ray 
analysis. We show the data in Fig. [8] along with three mod- 
els, each with M = 7 x 10"'' Mq yr~^ and l3 = 0.7. The best 
model has a constant /v = 0.08 above a radius Rd = 1.05 
Rt, below which the wind is assumed to be smooth. Note 
that this Rci is not necessarily the same as the X-ray onset 
radius, Ro. In fact, this small clumping radius, R^i is nec- 
essary to fit the data, as comparison with the model that 
assumes Rd = 1.3 i?, in Fig. |S] shows. In that model, the 
simulated strength of the core is much too low, reflecting 
the reduced Hq opacity (which scales as M^/fv) in the lower 



Figure 8. The Ha emission profile (black, solid) is well repro- 
duced by a model that includes clumping beginning at a radius 
of Rd = 1.05 Rt (red, dashed). The model with R^ = 1.3 Rt, 
(green, dash-dot) does not have enough emission in the line core. 
For comparison, we show a model with no clumping (blue, dot- 
ted), in which the model fails to produce any emission. All models 
have M = 7 X 10"'' Mq yr~^ and /3 = 0.7. Note that the narrow 
emission peak at line center likely has a significant contribution 
from nebular emission, which we do not model. 



wind. This early onset of wind clu mping has been found for a 
number of other stars as well tBouret et al ] |2005l : IPuls et al.l 
l200d) . We also note that the simplest case of a spatially con- 
stant /v actually reproduces the Ha data reasonably well for 
this star, in contrast to, e.g., A Cep, for which a rather stro ng 
radial dependence of /v is needed l|Sundqvist et al.ll201ll ). 



3.3 Global thermal modeling 

While the individual line profiles provide information about 
the kinematics of the shock-heated plasma, its spatial dis- 
tribution, and the degree of attenuation by the bulk wind 
in which the shock-heated plasma is embedded, complemen- 
tary information is contained in the overall spectral energy 
distribution. The broadband spectrum provides informa- 
tion about the temperature distribution of the shock-heated 
plasma and also the wind mass-loss rate via the degree of 
attenuation of the X-rays. So, to complement the spectral 
modeling of individual emission lines, described in the pre- 
vious two subsections, we fit global thermal models to the 
low-resolution zeroth-order spectrum of component A shown 
in Fig. H 

The model assumes optically thin, coUisional-radiative 
equilibrium ("coronal") emission {vapec, the (V ariable abun- 
dance) Astrophysical Plasma Emission Code l|Smith et al.l 
,200lh ). It is attenuated by the cold, neut ral interstel- 
lar m edium (via the xspec model tbabs l|Wilms et al.l 
I2OOOI )) and - for the EWS emission - the partially ion- 
ized stellar wind (via th e xspec custom model wmdtab^ 
(|Leutenegger et al] |2010| )V We model the EWS emission 
with a single isothermal vapec spectrum. To account for 



This model is also described on the windprofile custom model 
page on the XSPEC site. 
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the contribution of harder X-rays from the wind-wind in- 
teraction we include a second thermal emission componenlQ 
attenuated only by the ISM. 

This composite model is invoked in xspec as 
(vapecx windtabs + vapec)xtbabs. The free parameters of 
the model include the temperatures of the two vapec com- 
ponents, their emission measures, the characteristic mass 
column, E, = t^/k (g cm~^), of the wind absorption model, 
windtabs, and the interstellar column density. Fixed param- 
eters include the metallicity (fixed at solar, except for CNO) 
of the emission model and the wind velocity profile (de- 
scribed by /3 — 0.7) of the wind attenuation model . We 
note that the windtabs model (|Leutenegger et ai]|2010l ) has 
two features that make it distinct from interstellar attenua- 
tion models and make it more appropriate for the modeling 
attenuation by a stellar wind with embedded shocks: (1) 
it incorporates atomic cross sections from partially ionized 
species (e.g. O iv rather than neutral O) and assumes that 
H and He are fully ionized; and (2) it uses an exact ra- 
diation transport model appropriate to a n emitter spatially 
distrib uted within the absorbing medium. ^Leuteneg ger et al] 
show that the attenuation from this realistic wind 
transport model differs significantly from the exponential 
("slab") attenuation implemented in ISM absorption mod- 
els (their figures 5 and 8). For the fitting we r eport on here, 
we use a solar abundance (|Asplund et al.|[2009l ) wind opacity 
m odel, but with altered C, N, and O abundances according 
to lTaresch et al.l (|l997l ) , just as we did for the analysis of the 
ensemble of r* values discussed in the previous subsection. 
The vapec emission model assumes the same abundances as 
the windtabs absorption model. 

We fit the above-described composite thermal emission 
with wind-plus-ISM absorption model to the zeroth-order 
spectrum extracted from the seven coadded pointings. This 
low-resolution CCD spectrum has significantly better signal- 
to-noise than the dispersed grating spectra, but at a resolv- 
ing power (E/AE) of only a few tens. It extends to lower 
energies (0.5 keV) than the dispersed spectrum effectively 
does. The zeroth-order spectrum does not suffer from signif- 
icant pile-up effects (unlike bare ACIS observations of the 
same star). We use the goodness of fit statistic and 

put confidence limits on the fitted mod el parameters using 
the Ax^ formalism of lPress et all l|2007l ). 

We let the interstellar column density be a free param- 
eter of the fit, and found a best-fit value of A^h = 4.1 x 10^^ 
cm~^, which is very close to th e value implied by the color 
excess (Ah = 3.7 X 10^1 cm-2 ||Gagne et al.ll201ll )). The 68 
percent confidence limits for Ah (ISM) extend from 3.5 x 10^^ 
cm~^ to 5.0 X lO'^^ cm"'^. The best-fit model has emis- 
sion component temperatures of kT = 0.61 ± .02 keV and 
kr = 3.3llo'g9 keV, where the hotter component has only 
6 percent of the total emission measure. We conjecture that 
this hotter component represents a small amount of CWS X- 
ray emission, presumably associated with the non- thermal 
radio emission detected in the system jBenaglia et al,,2006 ') . 
Its contribution is negligible below photon energies of 2.5 
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Figure 9. The same zeroth order ACIS CCD spectrum shown 
in Fig. [2] here fit with a two-temperature apec thermal emission 
model (red histogram), where one temperature component (0.6 
keV) is attenuated by the stellar wind as well as the interstellar 
medium and the other (3.3 keV) is attenuated only by the ISM. 
Note the presence of strong Si xin emission just below 2 keV. The 
vast majority of the emission in this spectrum is line emission, but 
due to the low resolution of the detector as well as the presence of 
many weak, blended lines, the spectrum looks relatively smooth. 
The inset figure shows the same data with a model identical to the 
best-fit model, except that the wind absorption (S« in windtabs) 
is zeroed out. This model spectrum makes the significance of the 
wind absorption effect quite obvious. Nearly 80% of the emitted 
EWS X-rays are absorbed before they can escape from the wind. 



keV. This provides further confirmation that the wind-wind 
X-rays do not affect the line profiles we discussed in the 
previous subsection. 

The characteristic wind mass column density, E* in the 
windtabs model, is found to be E, = 0.0522lo;o^4g g cm"^. 
Using the wind term inal velocity and stellar radius from 
iRepolust et al.l l|2004l ). this corresponds to a mass- loss rate 
of M = 5.18t};4g X 10" Mq yr"^ The biggest contribution 
to the uncertainty on the wind mass column density is the 
uncertainty on the interstellar absorptiorjfl. The quality of 
the fit is good, with a reduced of 1.01 for 113 degrees of 
freedom. It is shown in Fig. |9] 

We note the consistency of the broadband mass-loss rate 
determination with the independent determination from the 
ensemble of line profile shapes, discussed in the previous sub- 
section. The wind attenuation is significant, as we show in 
Fig. [9l where we also include a model without wind attenua- 
tion. The X-ray luminosity (corrected for ISM attenuation) 
of the best-fit model is Lx — 7.29 x 10^^ ergs s~^ (giving 
Lx/LboI = 1-3 X lO"'^), but when we also correct for wind 
attenuation, we find Lx ~ 3.30 x 10'^'' ergs s~^. This im- 
plies that 78 percent of the X-ray emission above 0.5 keV 
produced by embedded wind shocks is absorbed before it 
escapes the wind (which can be seen graphically in the inset 
of Fig. O. Note that nearly all of the X-ray emission below 
0.5 keV will also be attenuated. 

To further test the plasma emission temperature result. 



* Although the temperature distribution of the shocked plasma 
is certainly more complex than two discrete temperatures, we find 
that adding more temperature components does not improve the 
fit quality. 



° When we fix the IS M column density at A^h = 2.5 X 10^^ cm"^ 
llTaresch et al.lll997^ . the wind mass column density, and along 
with it, the mass-loss rate, increases by a factor of 1.8. 



10 D. Cohen et al. 



we can examine the temperature-sensitive iron L-shell line 
complexes in the 10 to 17 A wavelength range of the grat- 
ing spectrometer. The data do not have good signal-to-noise 
in that wavelength region and likely also sufTer from modest 
contamination by dispersed photons from HD 93129B, as we 
discussed in §2. The contamination was primarily a concern 
with respect to the emission line shapes. But a contamina- 
tion level of 20 or 25 % (based on Fig. [2]) will not significantly 
skew the spectrum over a broad wavelength range. Because 
we want to ignore the contaminated line shapes, and to en- 
hance the signal-to-noise ratio, we rebinned the MEG spec- 
trum, using a 20 counts per bin criterion. We then fit the 10 
to 16 A portiorjfl of the MEG spectrum with an apec emis- 
sion model, including wind attenuation via windtabs and 
interstellar absorption via tbabs. We fixed the wind mass 
column parameter of windtabs at the best-fit value found 
in our global fitting reported on earlier in this subsection 
and fixed the interstellar column density at A^h = 3.7 x 10^^ 
cm~^. 

Because of the large number of charge states with 
closely spaced ionization energies, Fe L-shell emission is 
a sen sitive diagnostic of plasma temperature (Bchar ct al.i 
I2OOII ). Ne-like Fe xvil, with strong emission lines near 15 
and 17 A, is present over a moderately wide range of tem- 
peratures, and dominates over higher charge states below 
about kr = 0.3 keV. Fe xvill has strong lines near 14.2 
and 16 A and Fe xx has strong lines near 12.8 A. Higher 
ionization states have strong lines primarily between 11 and 
12 A. 

In Fig. [To] we show the binned MEG data along with the 
best-fit thermal emission model, which has a temperature 
of kr = 0.58 keV, in very good agreement with the zeroth 
order spectral fitting. Here, though, we can see which specific 
lines are and are not contributing to the observed flux. The 
strong emission of Fe xvil is seen clearly near 15 A, but there 
is no strong emission from the higher charge states. Thus, 
the dominant plasma temperature is constrained to be about 
0.6 keV or lower, while only very small contributions from 
higher temperatures are compatible with the data. 



3.4 HD 93129B 

The zeroth order spectrum of HD 93129B, the 03.5 com- 
panion at a separation of 2.7", is shown in Fig. [2] We have 
also flt it with a thermal emission model with both wind 
and ISM absorption (but no CWS component), much as 
we did for HD 93129A. Because of the low signal-to-noise 
of this spectrum, we held the interstellar column density 
fixed at A^H = 3.7 x 10^^ cm~^, corresponding to the mea- 
sured color excess. Our best-fit model has a temperature of 
kr = 0.31 ± .06 keV and a wind absorption mass column of 
E, = 0.0221 j^ Qjjg g cm~^, corresponding to a mass-loss rate 
of M « 1 X 10~® Mq yr~^, with a factor of two uncertainty. 

We note that this source emi ssion tempe r ature is some- 
what lower than that reported in lNaze et al.l (|201ll ). This is 
most likely due to our inclusion of wind absorption, which 
hardens the emergent spectrum. The ISM-corrected X-ray 
luminosity of our best-fit model is Lx = 1.5 x 10^^ ergs s~^. 



^ Absorption renders the iron line complex near 17 A undetected. 
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Figure 10. The binned MEG spectrum (squares are -1 order 
and circles +1 order) is shown along with the best-fit (MEG -1 
order) model (red). Cooler (0.3 keV, green) and hotter (1.5 keV, 
blue) models are also shown, making it clear that the presence 
of Fe XVII emission near 15 A in the data and the low levels of 
emission near 11 A from higher ionization stages are incompatible 
with very hot (kT > 1 keV) plasma from colliding wind shocks. 
Ne IX and X emission complexes are also visible at 13.5 A and 
12.1 A, respectively. Their ratio is temperature dependent and 
consistent with temperatures between 0.3 and 0.6 keV. We indi- 
cate the wavelengths of some of the strongest lines at low and 
at high temperatures. Note that there are other, low-temperature 
lines of, e.g., Ne IX contributing to the observed flux in the short 
wavelength side of the spectral region displayed in this figure. 



corresponding to Lx/LboI = 1.2 x 10 
value we find for component A. 
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very similar to the 



4 DISCUSSION 

The Chandra observations of HD 93129A are consistent 
with the embedded wind shock scenario that is generally 
assumed to apply to all O stars. The X-ray luminosity is 
Lx = 1.3 X 10~^Lboi; very much in line with the canonica l 
value for O star wind-shock sources ijPallavicini et al1ll98ll ) , 
although it is roughly a factor of two higher than the aver- 
age value found for O stars in Carina (|Naze et al ][20Tlh . The 
emission temperature is quite low (0.6 keV), as is expected 
from the EWS mechanism, although there is evidence for 
a small ('^ 6%) contribution from hotter plasma (3.3 keV). 
The presence of this hotter component is not surprising, 
given the detection of non-thermal radio emission from the 
system, presumably associated with colliding wind shocks 
involving components Aa and Ab, separated by roughly 100 
AU. But we stress that the colliding- wind contribution to 
the overall X-ray spectrum is minimal. 

Although the EWS emission temperature of 0.6 keV is 
low, the overall spectrum is relatively hard. We have shown 
here, using both the individual line profiles and also the 
broadband, low-resolution CCD spectrum, that this is due 
to attenuation by the star's dense stellar wind. The two dif- 
ferent, largely independent, manifestations of the wind at- 
tenuation lead to consistent mass-loss rate determinations 
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of between 4.7 and 7.0 x 10^^ Mq yr~^, which represent a 
factor of several reduction in the mass-loss rate over tradi- 
tional values determined from density- squared diagnostics. 
As both X-ray mass-loss rate diagnostics are insensitive to 
density squared effects, and because some clumping is cer- 
tainly expected in the wind of HD 93129A (|Lepine fc Moffati 
l2008l ). the new, lower mass-loss rate seems quite reasonable, 
and is in line with the factor of three mass-loss rate reduc- 
tion seen in C Pup (|Puls et al.ll2006l : ICohen et al.lbOlOl '). This 
is verified by our modeling of the Ha line, which shows that 
M = 7 X 10~® Mq yr~^ provides a good fit if a constant 
clump volume filling factor /v = 0.08 is assumed. To ex- 
plain the strong core of the Ela line, a clump onset radius 
of Rci = 1.05 R* also must be assumed. We note that the 
X-rays are not produced at radii this small. Presumably the 
wind shocks associated with clumps at the slow-moving base 
of the wind are not strong enough to produce X-rays. 

The kinematics of the X-ray emitting plasma in IfD 
93129A determined from the line widths are consistent with 
the terminal velocity of the bulk wind {voo ~ 3200 km 
s~^), assuming a constant X-ray filling factor above some 
onset radius, Ro. We derive Ro values for each line or line 
complex from the profile fitting, and find onset radii that 
are somewhat lower than, but still statistically consistent 
with, the typical value of Ro ~ 1.5 R* from numerical sim- 
ulations of the line-driving instability (Owocki et al.l 1 19881 : 
iFeldmeier et al.lll997l : [r unacres fc Qwockill2002l ). However, 
if the wind terminal velocity is lower than we assumed for 
our fitting, then the Ro values would increase by 0.1 or 0.2, 
to _Ro ~ 1.5 R* (for an assumed Voo = 3000 km s~^). 

The location of the X-ray plasma is further con- 
strained by the forbidden-to-intercombination line ratios in 
the helium-like line complexes observed in the grating spec- 
tra. Using a model that includes both the broadening and 
attenuation effects on the line profiles in conjunction with 
the altered forbidden-to-intercombination line ratios due to 
UV photoexcitation, we find that all three helium-like line 
complexes are completely consistent with the X-ray emit- 
ting plasma being distributed throughout the wind of HD 
93129A, starting at a height of only several tenths of a stellar 
radius above the photosphere. Using a simpler model where 
we fit the forbidden-to-intercombination line ratio directly, 
using Gaussian profiles, we find low f /i ratios for the Si xill 
and the Mg xi complexes. For both these complexes, the ra- 
tio expected if photoexcitation is unimportant, as would be 
the case in a binary wind-wind collision zone far from ei- 
ther star's photosphere, is ruled out with greater than 95% 
confidence. 

Finally, we must discuss alternative interpretations of 
the spectra and, especially, the spectral lines. The profiles, 
while skewed and blue-shifted, are not as asymmetric as ex- 
pected given the very high mass-loss rate traditionally found 
for the star. Our interpretation is that this is due to an 
actual mass-loss rate that is modestly lower than the tradi- 
tional value, but in principle, it could also be due to porosity 
associated with optically thick and presumably large-scale 
wind clumping. However, not only are the required poros- 
ity lengths very high, but the same clumps that would need 
to be invoked to generate a porosity effect would also lead 
to a mass-loss rate reduction due to their effect on density- 
squared diagnostics. If, for example, the required porosity 
lengths were achieved with clumps having a size scale of 0.25 



R, and a uniform filling factor of fv =0.1 in the context 
of an Ha mass-loss ra te of M = 1.8 x 10~ ^ Mq yr"'^ that 
assumes no clumping (|Taresch et al.|[l997l ). then the filling 
factor alone, via its effect on the density-squared mass-loss 
rate, would reduce the mass-loss rate inferred from the Ha 
to a value consistent with what we find from the X-rays 
without having to invoke any porosity - as we have shown 
at the end of §3.2 and in Fig. [S] Any additional effect from 
porosity would make the mass-loss rates too low and the 
X-ray line profiles too symmetric. 



5 CONCLUSIONS 

We have shown that the Chandra grating spectrum of the 
extreme O star, HD 93129A, can be understood using the 
same paradigm that explains the canonical embedded wind 
shock source, ( Pup, once adjustments are made for its larger 
wind terminal velocity and mass-loss rate. Specifically, the 
kinematics of the X-ray emitting plasma are consistent with 
shocks embedded in a [3 — 0.7, Voo = 3200 km s~^ wind 
starting at several tenths R,, and that the attenuation sig- 
natures in the line profiles are consistent with a mass-loss 
rate of 4.7 to 7.0 x 10~^ Mq yr~^, representing a modest re- 
duction compared to traditional mass-loss rates determined 
from Ha measurements that ignore the effects of clumping, 
and showing consistency with Ha modeling that includes 
modest c lumping, in line with w hat is seen in LDl sim- 
ulations (|Dessart fc Owockill2005l ). This mass-loss rate re- 
duction of a fact or of three to four is consistent with that 
found for ( Pup l|Cohen et al.ll2010^ . We have also demon- 
strated for the first time that modeling wind absorption of 
X-rays for both line profiles and for the broadband spec- 
tral energy distribution leads to consistent results when a 
physi cally realistic model of the broadband wind attenua- 
tion (|Leutenegger et aLll2010l ) is used. 

The global spectral modeling indicates that the domi- 
nant thermal emission component has quite a modest tem- 
perature, of roughly 0.6 keV, as predicted by EWS mod- 
els. The observed overall hardness of the spectrum is at- 
tributable to wind attenuation, rather than high plasma 
temperatures. The low dominant plasma temperature is also 
manifest in the low Si xiv/Si xill ratio, which is consistent 
with the value found in the Chandra grating spectrum of 
^ Pup. There is likely a small amount of hard X-ray emis- 
sion from colliding wind binary interaction between com- 
ponents Aa and Ab, many tens of AU from either star's 
photosphere. Because of the large separation of the compo- 
nents, this X-ray emission makes a small contribution to the 
overall X-ray spectral properties, representing less than 10% 
of the system's X-ray luminosity. The helium- like f /i ratios 
also provide evidence that the bulk of the X-rays arise in 
embedded wind shocks. 

As HD 93129A is the earliest O star known, and has one 
of the strongest winds of any O star, the work presented here 
strongly suggests that the embedded wind shock scenario, 
as described by numerical simulations of the line-driving in- 
stability, is widely applicable to O stars, even those with ex- 
tremely strong winds. And that X-ray line profile analysis, 
especially in conjunction with broadband spectral modeling, 
provides a good means of making a dumping-independent 
mass-loss rate determination for O stars with dense winds. 
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